Nonlinear magnetoresistance of an irradiated two-dimensional electron system 
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Nonlinear magnetotransport of a microwave-irradiated high mobility two-dimensional electron 
system under a finite direct current excitation is analyzed using a dc-controlled scheme with photon- 
assisted transition mechanism. The predicted amplitudes, extrema and nodes of the oscillatory 
differential resistance versus the magnetic field and the current density, are in excellent agreement 
with the recent experimental observation [Hatke et al. Phys. Rev. B 77, 20f304(R) (2008)]. 
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The prediction 1 - and detection^ of radiation in- 
duced magnetoresistance oscillation (RIMO) in two- 
dimensional (2D) electron systems (ES), especially the 
discovery of the zero-resistance stated, have stimulated 
intensive experimental 6 ^^^^^^^^ and 
theoretica^^^^^^^^ 2 ^^^ 2 .^ 3 .^ 5 . studies 
on this extraordinary transport phenomenon of electrons 
in very high Landau levels. 

Despite the fact that basic features of RIMO have 
been established and the understanding that it stems 
from impurity scattering has been reached, so far 
there has been no common agreement as to the accu- 
rate microscopic origin of these giant resistance oscil- 
lations. Presented in different forms, many theoreti- 
cal models^^^^^ 2 ^^^^ 7 . consider RIMO to 
arise from electron transitions between different Landau 
states due to impurity scattering accompanied by ab- 
sorbing and emitting microwave photons. This origin 
is called the "photon-assisted transition" or "displace- 
ment" mechanism. A different origin, called the "inelas- 
tic" or "distribution function" mechanism,— ! 30 ' 31 consid- 
ers RIMO to arise from a microwave-induced nonequi- 
librium oscillation of the time-averaged isotropic elec- 
tron distribution function in the density-of-states (DOS) 
modulated system. Both mechanisms exist in a real 2D 
semiconductor and have been shown to produce magen- 
toresistance oscillations qualitatively having the observed 
period, phase and magnetic field damping. The "dis- 
placement" mechanism predicts a well-defined photore- 
sistivity with given impurity scattering and Landau-level 
broadening, while the "inelastic" mechanism yields an 
additional factor proportional to the ratio of the inelas- 
tic scattering time r- ln to the impurity-induced quantum 
scattering time r q M- Since the inelastic scattering time 
Tj n or the thcrmalization time 7th being the property 
of a nonequilibrium state and contributed by the direct 
Coulomb interactions between electrons and by all other 
possible impurity- and phonon-scattering mediated effec- 
tive electron-electron scatterings^ is very hard to de- 
termine theoretically or to measure experimentally, the 
sharp controversy whether T m /T q 3> 1 or Ti n /r g -C 1, 
i.e. which mechanism plays the dominant role in the 
experimental systems ) 23 ' 30 has been an unsolved issue. 
The detailed comparison between theoretical predictions 
and experiments may provide a useful way to distinguish 



them. 

Introducing additional parameters into microwave- 
illuminated 2DESs, such as dc excitations, can be of 
help to distinguish different models and mechanisms. It 
has been shown that a finite current alone, can also in- 
duce substantial magnetoresistance oscillation and ze- 
roresistance without microwave radiation i 38 ' 39 ' 40 ' 41 ' 42 ' 43 
Simultaneous application of a direct current and a mi- 
crowave radiation leads to very interesting and compli- 
cated oscillatory behavior of resistance and differential 
resistance! 36 ' 37 ' 44 Recent careful measurement s 45 ' 46 dis- 
closed further details of such nonlinear magnetotransport 
in a high-mobility 2D semiconductor under both ac and 
dc exitations, allowing a careful comparison with theo- 
retical predictions. 

Our examination is based on a current-controlled 
scheme of photon-assisted transport* 2 ^ which deals with a 
2DES of short thermalization time having N s electrons in 
a unit area of the x-y plane and subject to a uniform mag- 
netic field B = (0, 0, B) in the z direction. When an elec- 
tromagnetic wave with incident electric field E\ s sin cot 
irradiates perpendicularly on the plane together with a 
dc electric field Eq inside, the steady transport state of 
this 2DES is described by the electron drift velocity Vq 
and an electron temperature T c , satisfying the force and 
energy balance equations^ 3 - 

N s eE + N s e{v x B) + F = 0, (I) 
N s eE ■ v + S p - W = 0. (2) 

Here, the frictional force resisting electron drift motion, 

oo 

F = J2\U(q\\)\ 2 £ q [ |J 2 (e)n 2 (q || ,^ + ^), (3) 

qii n— — oo 

is given in terms of the electron density correlation func- 
tion LT2 (<7|| , SI) , the effective impurity potential U(q«), a 
radiation-related coupling parameter £ in the Bessel func- 
tion J n (C)i an( i = Q\\ ' v o- The electron energy absorp- 
tion from the radiation field, S p , and the electron energy 
dissipation to the lattice, W, are given in Ref.[23j. The 
nonlinear longitudinal resistivity and differential resistiv- 
ity in the presence of a radiation field are obtained from 
Eq. JTJ) by taking Vq and the current density J = N s evo 
in the x direction, Vq = (vq, 0, 0) and J = (J, 0, 0), 



-F /(Nyvo), r, 



~(dF /dv )/(Ny). (4) 
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FIG. 1: (Color online) Differential magnetoresistivity r„ vs 
e w = oj/ljc under fixed bias current densities from J — to 
J = 0.16 A/m, in 0.01 A/m increments. 
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FIG. 2: (Color online) (a) Differential resistivity r xx vs ej at 
fixed from 2 to 3.5 in steps of 0.0625. (b) Amplitude of 
r„ oscillations, (b) Positions of r xx maxima extracted. 



We have calculated the differential resistivity r xx from 
above equations (taking up to three-photon processes) 
under different magnetic fields B and bias drift velocities 
vq for a GaAs-based heterosystem with carrier density 
iV s = 3.7 x 10 15 /m 2 and low-temperature linear mobility 
/io = 1200m 2 /Vs at lattice temperature T = 1.5 K, irra- 
diated by a linearly x-polarized microwave of frequency 
lo/2tt = 69GHz with incident amplitude E- ls = 3.6 V/cm. 
The elastic scatterings are assumed due to a mixture of 
short-range and background impurities, and the Landau- 
level broadening T is taken to be a Gaussian form with 
a broadening parameter a = 7^ 

Figure 1 presents the calculated r xx versus e u = u)/oj c 



(uj c = eB /mis the cyclotron frequency) at fixed bias drift 
velocities from 2vo/vp — (vf is the Fermi velocity) to 
2 x 10~ 3 in steps of 1.25 x 10~ 4 , corresponding to current 
densities J = to 0.16 A/m in steps of 0.01 A/m. The 
J = case exhibits typical RIMO with a sequence of 
resistance maxima (l + ,2 + ,3 + ,4 + and 5 + ) and negative 
values around the resistance minima 1 _ and 2~. With in- 
creasing J to 0.16 A/m, the maxima 2 + , 3 + and 4 + (min- 
ima 2~, 3 _ and 4~) evolve into minima (maxima) having 
seemingly little change in the B positions. Further, all 
the curves cross approximately at e w = 1.5,2,2.5,3,3.5 
and 4.5, indicating that J at this range has little effect on 
photoresistance there. These and other features of Fig. 1 
reproduce what was exactly observed in Refl46l. 

Figure 2(a) shows the calculated r xx versus ej = u)j/u> c 
(ujj = 2kpVo, k-p is the Fermi wave vector of the 2D elec- 
tron system) at fixed e w from 2 to 3.5 in steps of 0.0625. 
Traces are vertically offset in increments of 0.625 for 
clarity. The amplitudes Ar of r xx -ej oscillations are max- 
imized around e w ~ 2.2,2.8 and 3.2 and strongly sup- 
pressed at e w ~ 2.5 and 3.5, as shown in Fig. 2(b).— The 
positions of r xx maxima extracted are plotted in Fig. 2(c) 
as dots. All these are in excellent agreement with the ex- 
perimental results [Fig. 2(a), (b) and (c) of Ref lioT] . 

In the present current-controlled transport model the 
oscillations of R xx and r xx are referred to the behavior of 
function IT2 (<Z|| , ujq + uuj) in Eq. ([3]). The electron density 
correlation function 112(911 > ^) is essentially a multiplica- 
tion of two energy-fi shifted periodically modulated DOS 
functions of electrons in the magnetic field^ 3 - Its period- 
icity with changing frequency Q — > Q + uj c at low tem- 
peratures and high Landau-level occupations, determines 
the main periodical behavior of magnetoresistance. 23 The 
previous examination 36 focused on the node positions of 
the oscillatory peak- valley pairs of R XXl which appear 
periodically roughly along the lines e w + V e j = m — 
1, 2, 3, 4, ... in the e^-e^ plane, where e w = lo/lo c is the con- 
trol parameter of RIMOs, and €j = ujj/u> c is the control 
parameter of current-induced magnetoresistance oscilla- 
tions, and rj ^ 1, dependent on the scattering potential. 42 

The maxima of differential resistivity r xx show up at 
lower values in the €j axis in comparison with the node 
positions of related valley-peak pairs of R xx , and its ap- 
pearance exhibits a periodicity Aej < 1.— In the e^-ej 
plane, the differential resistance maxima are expected to 
show up roughly in the vicinity along the lines 



+ Ae, 



1,2,3,4, 



(5) 



with A > 1, dependent on the scattering potential (A = 
1.04 for the system on discussion). Eq. qualitatively 
accounts for the periodical change of r xx in a large scale 
in steps of A(e w + A ej) = 1. 

Under strong microwave irradiation, as in the present 
case, the role of virtual photon process [the n = term 
in the sum of Eq. ([3])] is negligible due to samll Jq (£),— 
and main contributions to resistivity come from n = 
±1, ±2, ... terms (single- and multiple-photon processes). 
Noticing that the frequency differentiate TE^ , S^) is an 



3 



even function of Q and considering contributions from 
scatterings parallel and antiparallel to the drift velocity 
Vq and from ±|n| terms, we see that, in the case of fi- 
nite bias current, the r xx behavior is determined by the 
sum of two terms: (a) n 2 (<7||, \n\uj + q^v cos8) and (b) 
n' 2 (q||, \n\uj — q\\vocos9). Depending on the n 2 (<7||,fi) 
function behavior in the vicinity of £1 = \n\ui, effects of 
these two terms can be cancelled or added, completely 
or partly, at different locations of LT 2 (<7||,f2) func- 
tion reaches maxima (positive) at around H/oj c = N — j, 
reaches minima (negative) at around H/oj c = N + j, and 
passes through zero (changing sign) at around Q/u c = N 
and N + \ for all integers N > 2. Thus, at e w ~ I or 
1+2 (I = 2, 3, 4, ...) with which all involved \ti\uj frequen- 
cies are located around Nlu c or (N + ^)u> c , contributions 
from (a) and (b) are almost cancelled out for modest Vq. 
In the case of e w ~ I — \ [e u 1+ |], there always exists 



a term of frequency (N — j)oj c [{N + \)u> c ] in \n\ui, and 
contributions from (a) and (b) are positively [negatively] 
additive. These clearly account for the suppression of the 
current effect at e w ~ I and I + |, and the enhancement 
of it around e w ~ I — j and e w ~ I + j . 

Above discussions are general. The accurate behavior 
of resistivity r xx inside a period scale is relevant to the 
detailed shape of the DOS function. Figs. 1 and 2 repre- 
sent the result of a Gaussian-type DOS. The good quan- 
titative agreement with experiment without adjusting 
parameters indicates that the present current-controlled 
scheme of photon-assisted transport captures the main 
physics of RIMOs in the discussed quasi-2D system. 
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